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Abstract 

An analysis procedure is proposed to measure the antiquark distributions 
in the proton over the region 0.01 < x < 0.1. The procedure involves the 
measurement of high pt asymmetric direct photon and jet final states in pp 
interactions. This measurement can be made at the RHIC collider running 
in pp mode at an energy of -y/i = 500 GeV/c. This analysis identifies a re- 
gion of phase space where the contribution from quark-antiquark annihilation 
uncharacteristically approaches the magnitude of the contribution from the 
leading process, quark-gluon Compton scattering. The forward-backward an- 
gular asymmetry in the parton center of mass is sensitive to the antiquark 
content of the proton and the u parton density function can be extracted. 
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I. INTRODUCTION 



The distribution functions of quarks and gluons in a proton play a crucial role in the self 
consistent description of all hard inclusive proton interactions. The cross sections in both 
lepton-proton and proton-proton high-^ interactions are derived from the fundamental 
perturbative QCD parton interaction cross sections, folded with the parton distribution 
functions and integrated over available phase space. The main sources of experimental input 
for the parton distribution functions (or structure functions) include deep inelastic lepton- 
nucleon scattering and hard hadronic inclusive interactions with final states including high 
mass leptons pairs, high-p t direct photons, high-p t jets, W, Z or heavy flavors. 

Knowledge of the quark and gluon distributions of the proton has greatly increased in 
recent years. The first input has come from the high quality data on deep-inelastic structure 
functions in lepton-nucleon scattering [|],@. These experiments are generally most sensitive 
to valence quark distributions. They probe electromagnetically and are only directly sensi- 
tive to charged partons rather than gluons. The deep-inelastic scattering experiments have 
been less conclusive in the determination of the antiquark distributions. Analyses of data 
to extract the antiquark distributions usually involve various assumptions ||. 

Progress in determining the sea quark and gluon distributions in the proton have con- 
centrated on fits to a variety of types of data. The CTEQ Collaboration ||-^| and MRS 
group [|7] have made global fits to available experimental data. But, the determination of the 
antiquark distributions with hadronic interactions using fit procedures has been difficult be- 
cause most high-pf hadronic reactions are either: dominated by contributions from valence 
partons and gluons; involve too many contributing subprocesses to analyze unambiguously; 
or have a very small cross section in the high-pt region. Furthermore, much of the recent 
high energy collider data involves pp interactions where the most probable interactions tend 
to involve gluons and valence quarks rather than non-valence quarks. 

It is unclear how well the antiquark distributions have been determined in the fits by 
the CTEQ and MRS groups. Even when considering similar data sets, these groups find 
that their fits differ in the low x region ||J]. For example, at Q 2 = 25GeV 2 the u(x, Q 2 ) 
distribution from MRS at x ~ (0.01 — 0.02) is higher by (10 — 15)% than the same distribu- 
tion of the CTEQ Collaboration. For the d antiquark d(x, Q 2 ) distribution, this difference 
increases up to (15 — 20)% at same region of x and Q 2 . It has been difficult to pin down the 
errors associated with these fits. 

To verify and improve the results of these or other fits, measurements which are more 
directly dependent upon and more sensitive to the antiquark distributions are needed. In 
this paper, one such experiment and analysis procedure is proposed. A kinematic region is 
defined for high-p 4 direct photon and jet production from pp interactions in which the uu 
annihilation process is uncharacteristically important. The region involves parton pairs with 
asymmetric momentum fractions and the analysis concentrates on the angular distribution 
of the parton-parton interaction. 

In general, the dominant subprocess for high-p t direct photon and jet production is the 
u quark-gluon Compton process. Since the direct photon cross section is proportional to the 
square of the quark or antiquark charge, the d and s quark contributions are suppressed by a 
factor of four relative to that of the u quarks. Thus, u-\-G^u + , yis the primary subprocess 
and the next most important process is u and u quark annihilation, u + u — > G + 7. The 
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u quark contribution to direct photon and jet production is enhanced over the d or s quark 
contributions both because of this factor of four and because of the reduced number of d or 
s quarks in the proton relative to u quarks. 

While the annihilation parton cross section is comparable in magnitude to that of the 
Compton process, the ratio of the distribution functions u{x) to G(x) is typically 1:10 in the 
region of x ~ 0.01—0.1. Generally, the u and u annihilation contribution to the direct photon 
cross section is only about 10%. These two types of parton cross sections do, however, have 
a different angular distribution. 

In the parton model, the direct photon-jet cross section is the sum over subprocess 
terms, each term is the product of a parton level cross section and the two associated parton 
density functions. The lowest order partonic cross sections for annihilation and Compton 
interactions are 

da 8 (i 2 + u 2 \ 

da 1 fs 2 + u 2 \ , . 

where s, t and u are the Mandelstam variables for the parton scattering subprocess. 

These two cross sections have very different angular distributions. The Compton cross 
section becomes very large for scattering angles corresponding to the u quark scattering 
backward relative to the incident u quark direction, sending the photon forward. The 
annihilation cross section becomes large for the photon emerging in the forward or backward 
direction. The idea of this analysis procedure is to focus on 7 + jet production in the region 
of the parton center of mass angular distribution, where the annihilation process is enhanced 
because of the small t singularity in the low order scattering graph. 

The variables X\ and X2 are defined to be the momentum fractions for the "first" and 
"second" protons respectively, and f? 7 as the angle between the photon and the "first" proton 
momentum direction in the center of mass of a partonic hard scattering subprocess. Exact 
definition of "first" and "second" protons will be explained in section |m[ From observable 
rapidities r^ 7 , i]j et , the transverse momentum p t and the pp center of mass energy s/s, the 
quantities Xi, x 2 and # 7 are defined as: 

Xl = 2L( e v^ +e iM) (3) 

x 2 = -^=(e~^ + e-"'-) (4) 
V s 

- exp (VM^h.) . (5 ) 




The forward-backward asymmetry A(xi,X2) is then defined as the ratio of the forward 
to backward cross section: 
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(6) 



The choice of cos($) = ±0.7 as the boundary defining the forward or backward regions, is 
somewhat arbitrary. 

For the two main subprocesses, annihilation and Compton scattering, the forward- 
backward asymmetry will have a different magnitude. The pair annihilation subprocess 
Aqq(xi,X2) = 1, and the Compton scattering case A qG (xi,x 2 ) ^> 1. Using this fundamental 
difference between A qq and A qG , the relative strengths of the q{x) and G(x) distributions in 
the proton can be determined. 



II. ANALYSIS OF ANGULAR ASYMMETRY 

This analysis considers the direct photon and jet production cross section from pp colli- 
sions in a small kinematic region defined by x± ± Ax and X2 ± Ax where x\ > x 2 . The asym- 
metric condition x\ > x 2 is important, enabling the observation of the forward-backward 
asymmetry of the direct photon and jet angular distribution. The proton antiquark density 
u(x 2 ) will be extracted from the direct photon forward-backward asymmetry in pp collisions 
as introduced in equation || 

For illustration, consider the region x\ = 0.175 ± 0.025 and X2 = 0.075 ± 0.025. The 
cos($ 7 ) distributions for the two main subprosses from a PYTHIA simulation of the direct 
photon and jet production in pp interactions [10] are plotted in figure [j]. The pp center of 



mass energy is ^fs = 500 GeV and the events are selected to be in the kinematic region 
defined by p t > 10 GeV, — 1 < r/ 7 < 2 and —1 < r]j et < 2 for both the jet and the photon. 
In the region cos($ 7 ) > 0.7 the contribution from the Compton scattering subprocess qG 
is more then five times greater than the contribution from the pair annihilation subprocess 
qq. In contrast, for cos(i? 7 ) < —0.7 the contribution from the gG-subprocess is only about 
twice the gg-subprocess contribution. 

The two dimensional plots of the event distributions under the conditions cos(i9 7 ) > 0.7 
and cos($ 7 ) < —0.7 are plotted in figure § for the same region of x\ = 0.175 ± 0.025 and 
x 2 = 0.075 ±0.025. From the results presented in figure ^|, we see that the cos(i9 7 ) > 0.7 cut 
selects large photon rapidity and small jet rapidity while the cos($ 7 ) < —0.7 cut selects the 
reverse ordering. Refining the definition of the asymmetry variable to include experimental 
acceptance, A(x%, x 2 ) is to be the accepted cross section (or events number) for cos(i9 7 ) > 0.7 
divided by the accepted cross section cos($ 7 ) < —0.7 for a particular (xi,x 2 ) bin. 

In terms of the parton distribution functions i(x) for parton % and the subprocess dif- 
ferential cross sections da it j / d cos(i? 7 ) for the process i + j — > jet + 7 , the asymmetry 
is 

T ! i(x 1 )j(x 2 )a^ j EN^ j (x 1 ,x 2 ) 

A(Xl ' X2) = j(x 2 ) a". = N~ j (x 1 , x 2 ) ' (7) 

i,3 hi 

where N^(xi,x 2 ) = i{x\) j(x 2 ) uf^ and i, j = u,u,d,d, s,s,c,c and G. The cross-sections 
af A are defined as: 
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cos(tf 7 )>0.7 

da. 



dcos(i?-v) 

cos(i? 7 )<-0.7 



dcos(0 7 ). (9) 



In an actual experiment, the divergence at the ±1 limits of integration will be removed 
by acceptance and experimental cuts. 

It is possible to estimate the asymmetry in the limit of full acceptance, with cross sections 
dominated by the small t and small u singularities of equations |8] and |9|. One can consider 
only contributions from the u and u quarks. A more complete analysis would also include 
small contributions from d, d, s, s and c, c pairs. With these approximations, the asymmetry 
is given by 

A(x x s Kc + N L + K* + N tu + N+ G + 

where = N^j(xi,x 2 ), i,j — u,u and G. Introducing variables for the ratios of u and u 
quark distribution functions u(x) and u(x) to gluons distribution function G(x), 

u(x) u(x) 
r{x) = , R(x) 



G(x) ' v ' G(x) 

and for the ratio of the integrated Compton cross section to the annihilation cross section, 
separately for the forward (+) and backward (— ) angular regions, 

B ± = ^. (11) 

Noting that a+ u = a~ a and cr+ G = a Gu = <r£ G = cr Gu , the approximate expression for the 
asymmetry reduces to: 

A(x 1 x 2 r(x 2 )) = c 2 r ( x 2) + cir(xi) + c ^ 
d 2 r(x 2 ) + dir(x\) + d ' 



where the constants q and di, i — 2, 1, and 0, defined in equation 12, are determined from 



the well known parton cross sections and u quark and gluon structure functions: 

c 2 ~ R(xx) + B~, d 2 ~ R{x x ) + B + 

Cl ~R(x 2 ) + B + , d 1 ^R(x 2 ) + B- 

c ~ B + R(x 1 ) + B-R(x 2 ) , d ~ B~R(xx) + B + R(x 2 ) 

For full acceptance, with cross sections dominated by the very large cos($ 7 ) = ±1 upper 
limit of integration, the integrals defined in equations ^| and |9] yield B + = 3/8 and B~ = 0. 
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By restricting the upper limits of integration to cos(-$ 7 ) = ±0.8, a more reasonable limit for 
a measurement at RHIC, B + increases to about 0.45 and B~ increases to about 0.1. Using 
these higher values, it is possible to estimate the size and range of observable asymmetries. 

Choosing X\ to be in the region where valence quarks dominate the structure func- 
tion [x\ ~ 0.2; R(xi) ~ 0.96 ) and choosing 22 in the region where the gluon distri- 
bution dominates (22 — 0.05; R(x 2 ) — 0.25 ), the observable range of asymmetry is 
2.2 > 74(21,22) > 0.75. This corresponds to r(22) ranging from zero to infinity while 
holding r(x\) at a fixed nominal value. For a nominal value of r(22) = 0.088, ^(21,22) 
would be about 1.65. 

It should be a very general conclusion that the asymmetry will be a ratio of linear 



functions of r(x\) and r(22) as shown in equation The values of the constants q and di 
will be evaluated with the PYTHIA Monte Carlo program for various (21, 22) bins, taking 
into account the experimental acceptance and p t cutoffs (see section p| ). It is noted that 
the dependence of the asymmetry upon the detailed model of the acceptance appears to be 
modest. The ratio of asymmetries for full acceptance (0.7 < | cos($ 7 )| < 1.0) to that for the 
limited acceptance described above (0.7 < |cos(i? 7 )| < 0.8) is only about 3:2, even as the 
formal cross sections range from infinite to finite. The range of asymmetries for the cuts 
discussed in the beginning of this section extends from about 0.75 to 2.2 reflecting "14(22) 
densities ranging from infinity to zero respectively. 

It is an important feature of this analysis procedure that the asymmetry is relativitely 
insensitive to some of the problems which have made the analysis of direct photon data 
difficult. Accurate modeling of the detector for a rapidly rising cross section near the edge 
of acceptance in cos(t? 7 ) has been one such problem. In this analysis, it is only the ratio of 
forward to backward efficiencies that must be understood. The price paid for this feature is 
that the inputs and outputs are not the quark or antiquark structure functions directly but 
are the ratios of those structure functions to the gluon structure function. 

With input from the ratios of quark to gluon structure functions and the ratios of in- 
tegrated parton Compton cross sections to quark-antiquark annihilation cross sections, the 
constants Cj and di are determined for various (21, 22) regions. These, along with the mea- 
sured asymmetry ^(21,22), will result in a determination of the antiquark to glue ratio. 
More specifically, this will determine a well defined combination of r(22) and r(xi), with a 
rather small coefficient for r(x\): 

r(x 2 ) + er(xi) = r £ , (13) 

where e and r e are determined from the constants q, di and the value of asymmetry A(x\, 22): 

ci - d 1 A(x 1 ,x 2 ) c - d A(x 1 ,x 2 ) 

d 2 A(x 1 ,x 2 ) - c 2 ' d 2 A(xx, x 2 ) — c 2 

By considering different regions (21,22), each providing a linear equation in r(2i) and 
r(22), unambiguous analysis should be possible. Unless the ratio r(x) grows rapidly with 
2 in the region 22 < 2 < 21, in contradiction of the conventional structure function sets, 
the effect of the r(2i) term is small, e ~ 0.1 (see section p| ). It is in this sense, that the 
asymmetry A{x\ 1 x 2 ) is a rather direct measurement of r(22) = u(x 2 )/G(x 2 ). 
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III. PYTHIA SIMULATION 



In this section we present a discussion of the p + p — > 7 + jet process at y/s = 500 GeV 



based on the PYTHIA simulation code [10|. The kinematic cuts are pt > 10 GeV, — 1 < 



77 7 < 2 and — 1 < 7]j et < 2, where /7 7 and rjj e t are pseudo-rapidity of the direct photon and 
jet, respectively. The p t cut is high enough to provide a fairly clean measurement of a direct 
(hard) photon with a managable background from tt° and 77 meson decay. The 77 7 and r]j et 
pseudo-rapidity cuts are consistent with the STAR detector acceptances assuming a Barrel 
and an Endcap Electro-magnetic Calorimeter (EMC). 

With one Endcap EMC, the STAR detector acceptance in pseudo-rapidity has a natural 
asymmetry - the position of the Endcap EMC in the pseudo-rapidity axis 1 < 77 < 2. Thus, 
the symmetry between two proton beam directions is broken and the side with the Endcap 
is defined to be at positive displacement from the interaction point and to have acceptance 
for positive rapidity. The convention will be that X\ is defined to be a momentum fraction 
of partons in the beam which moves in the positive direction, consistent with equations 0-[5|. 

The cross section for the direct photon plus jet production in the kinematic region defined 
by these rapidity cuts are shown for various [x\,x%) bins in table |. In table [TI| and table 
I II] the cross sections from the two main subprocesses, pair annihilation and the Compton 
scattering, are shown. The total direct photon and jet production cross section over this 
region is 9.3 nb. 

Figure ^ shows the accepted cross section distribution of x\ and x 2 . As shown in figure 
|3|, the pair annihilation contribution to the direct photon and jet production cross section 
is less than 10% (see, also, tables |IJ and |ITlf ). Two dimensional distributions of x\ and x 2 
are shown in figure 4. The asymmetry between x\ and x 2 reflects the asymmetric pseudo- 
rapidity cuts: —1 < rij e t < 2 and — 1 < ?7 7 < 2 , reflecting the STAR Barrel + Endcap 
coverage. 

The four small squares in figure [| are the regions analyzed in this paper. The values 
of constants q, di normalized by d 2 = 1, e and values of asymmetry A(xi,x 2 ) from these 
four regions are presented in table fV[ Also presented in table [V] is the accuracy of the 
predicted asymmetry A(xi,x 2 ) from PYTHIA Monte Carlo simulation code. The cos(i? 7 ) 
distributions of the events from one of these four regions is presented in figure |l|. 

The minimum value of x 2 , consistent with the cuts mentioned above in rapidity, p t and 
cos(?9 7 ), is x 2 nm ~ 0.02. The choice of a cos(i9 7 ) limit value of 0.7 is a trade-off between the 
need for a large asymmetry A(x\, x 2 ) and for a sufficiently large cross section. Reduction of 
the pt limit would provide coverage in regions of x 2 < 0.02. Lower x 2 can also be reached 
by decreasing the cos($ 7 ) limit. If asymmetry was defined by |cos($ 7 )| > 0.5, there is 
sensitivity down to x^ %n ~ 0.01 with the same p t cut. Such a choice involves larger accepted 
cross sections but weaker dependence of asymmetry upon the anti-quark ratio. 

The dependence of the forward-backward asymmetry upon r(x 2 ) is presented in figure 
|5|. The point at r(x 2 ) ~ 0.088 corresponds to the nominal CTEQ2L structure function sets 
which were used in PYTHIA (see figure ^). For the four (xi,x 2 ) regions considered here, 
when the nominal r(x 2 ) ~ 0.088 is assumed, the calculated asymmetries range from about 
1.5 to 2.0. 
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IV. ANALYSIS OF ERRORS 



To compare event rates for a variety of experiments at RHIC involving proton-proton 
colliding beams, RHIC experimentalists have agreed to consider a standard set of integrated 
luminosities and energies. At center of mass energy y/s = 500 GeV, the standard luminosity 
of 800 pb^ 1 is used to evaluate the expected statistical error. The fractional statistical error, 
SA/A, for the four regions considered is presented in table |I^. The typical statistical error 
in measurement of SA/A will be about 1%. The corresponding error introduced in the 
determination of 5r(x2) /r(x2) will be in the (3-7)% range (see figure 0), depending upon 
the actual value of r{x2)- 

Systematic errors will be larger than statistical errors. The systematic errors in the 
determination of r(x2) will come from various sources. The important systematic error 



sources will involve uncertainties in the determination of the constants in equations [12] and 
|T3"| and in the interpretation of the various sources of background. The following sources of 
errors have been considered. 

1. Contribution from partons other that u quarks and gluons G. The effect on 
asymmetry from all the other quarks and antiquark processes is small and calculable. 
Using the standard CTEQ2L, the scale for the change in asymmetry A(xi,x 2 ) due 
to the inclusion of other quarks and antiquarks is about 5 %. Even approximate 
simulation of these contributions would introduce a small uncertainty, 5 A/ A < 1%. 

2. Uncertainty in determination of q and di. The inputs to this analysis are the 
ratio of up quark to gluon structure functions R(x). It will be assumed that the ratios 
can be determined to 10%. As seen in figure [7], this implies an error in the (15-20) % 
at X2 — 0.025 and (30-40) % at x<i — 0.075 range for r{x-2). As mentioned above, the 
sensitivity to the detailed acceptance model near the limits of acceptance in rapidity 
and p t is not great. 

3. Effects from beyond Leading Order calculations. The main effect of the higher 
order corrections will be to increase the cross section by a A"-factor. While this is a 
large increase in cross section, the asymmetry A(xi,X2) is insensitive to the normal- 
ization. A full higher order calculation would determine a higher order correction to 
A(x\,X2) and the analysis procedure can be modified to take this into account. 

It was seen that in a Next to Leading Order (NLO) calculation of the CDF direct 
photon cross section |Hf , the change in the cos^) dependence is only at the 10% 



level. This suggests that the sensitivity of A(xi, X2) to a full NLO calculation may be 
rather modest. 

Background from jet fragmentation to high— p t photons. The standard back- 
ground to direct photons is from unmatched photons from ir° decays. Techniques 
have been developed to correct for this background It will be assumed here that 



the cross section of this background source can be directly measured and subtracted. 
Techniques for this analysis are well developed and will not be discussed here as a 
source of error in this analysis. 
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The background from photons produced directly in the jet fragmentation process is 
more difficult to estimate and subtract. The CDF analysis of direct photon plus 
jet production in pp interactions shows that the angular dependence is relatively en- 
hanced, in the forward-backward angle regions, over the NLO prediction. The angular 
distribution is enhanced beyond that predicted in NLO by about ~ (15 — 30)% in the 
region cos($ 7 ) > 0.5 ljl3| . They attribute this discrepancy to an underestimate of 
the fragmentation function of jets into photons, a calculation done only in leading 
order. This is about the level of modification expected for a full NLO calculation of 
jet fragmentation. 

It is reasonable to assume, that the contamination of the direct photon signal from jet 
fragmentation will be symmetric in cosf^). This assumption can be checked by a full 
NLO calculation of jet fragmentation into photons. For example, if the contamination 
was about 30% and symmetric in cos^), A{x\, x 2 ) would be reduced from about 2.0 
to about 1.63. If the correction was known to be (30 ± 10)%, the error would be 
5 A/ A < 7%. 

5. Effects from k t . A source of theoretical and experimental discrepancy is the initial 
state transverse momentum prior to a hard parton collision. An analysis of direct 
photon data fl4| suggests that this k t distribution is broader than expected, with a 
width as large as (2 — 3)GeV/c in this kinematic region of interest. Folded with the 
rapidly falling p t dependence, this results in an excess cross section at the lower limit 
of pt- In an analysis which anticipates a cross section proportional to the structure 
function, this is a great complication. 

The crucial point here is that it is not u(x) which is measured in this analysis but 
the ratio r(x) = u(x)/G(x). This observable is likely to be a more slowly varying 
function of x than the structure functions themselves in the small x region. For the 
analysis described here, a 50% increase in the cross section due to these effects would 
not be much of a problem. The main effect would be to increase the uncertainty of X2, 
somewhat smearing the bin over which the ratio r(x 2 ) is determined. The resulting 
error in x 2 is 5x 2 /x 2 — (10 — 20)%. We have shown that the asymmetry A(xi,x 2 ) 
varies rather slowly with x 2 for the cases considered. It is noted in equation |3| that 
the ratio of x\/x 2 is well determined by the rapidity and is not very sensitive to k t 
smearing. 

It appears, that with reasonable expectations of an eventual complete set of NLO calcu- 
lations, the largest source of error will come from the uncertainty in R(x) = u(x)/G(x). It 
does not seem too optimistic to expect a determination of r(x) = u(x)/G(x) to an accuracy 
of (15 — 30)%, at several values of x in the 0.02 to 0.1 region, using this method of analy- 
sis. The uncertainty in this measurement of the u structure function will be limited by the 
uncertainty in the u quark and gluon distributions. 
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V. SPIN DEPENDENCE 



The greatest interest in pp measurements at RHIC involves measurements with polar- 
ized beams. With longitudinally polarized beams, the dependence of the cross section on 
beam polarization can be measured to determine the polarization of the quarks and glu- 
ons within the proton. Another asymmetry, All, is defined to be the difference between 
the cross section for spin aligned and anti-aligned protons. This quantity is sensitive to 
the polarized parton distributions. A region is considered which is centered about (xi,x 2 ), 
with an unpolarized cross section a(xi,x 2 ) given by the average of two polarized cross sec- 
tions cr^(xi,x 2 ) and a^(xi,x 2 ), corresponding to longitudinal spins aligned or anti-aligned 
respectively. All(xi,x 2 ) is: 

A ( rv> ry \ ^ ^ {X\ ,X 2 )~ ^ 1 ^ ,X 2 ) 

^-LL\X 1 ,X 2 ) — CT TT( a;i ,a; 2 )+ CT Ti(a: 1 ,a;2) 

(14) 

= Eij P l L (xi)P J L (x 2 )f iJ (x 1 , x 2 )afj (cos(^)) 

where fij(xi,x 2 ) is the ratio of the unpolarized cross section from partons % and j to the 
sum of all processes in the region (xi, x 2 ), P^x) is the longitudinal polarization of a parton 
of type i at momentum fraction x relative to the longitudinal spin of the proton, and afj is 
the parton level analyzing power for this subprocess which can depend on the parton center 
of mass scattering angle. For annihilation and Compton scattering, the analyzing power is 



0: 



<>% = 4 L = - 1 ( 15 ) 



and 



s 2 - u 2 



1 - cos(^ 7 )) 5 



%G ~ %G ~ ~ 2 + ^ ~ 4 + (1 _ cos( ^ ))2 ( 16 ) 



LL_ LL _ S 2 -t* _ 4 - (1 + CQS(^ 7 )) 2 

a Gq - °>G, ~ §2 ~ $ ~ 4 - (1 - cos( ^ ))2 [17) 

If the magnitude of the anti-up quark structure function is near the nominal value 
r(x 2 ) ~ 0.1, then it is interesting to consider A\ L and A~[ L defined from considering cases 
of cos(-# 7 ) > 0.7 and cos(-# 7 ) < —0.7 respectively. Considering only the two most important 
terms in the sum from equation |14], 

A L l(x 1 ,x 2 ) ~ P u (xi)P G (x 2 )f UiG (xi,x 2 )afy(cos('d 7 )) 

(18) 

+P u (Xi)Pu(x 2 )f Ut u(xi, X 2 )(-l) 

Many of the factors in equation |l8| are to some degree known. In the region of x\ considered 
here, P u {x\) — 0.3 [p^ , p^j . The functions fij(xi,x 2 ) are the fractional contribution from a 
process involving % and j in the initial state and can be evaluated in PYTHIA. 
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The result is that in the forward region, the gluon polarization which is also measured 
in other kinematic regions will be nearly proportional to A\ L . However, in the backward 
photon region, A~[ L will be as sensitive to the anti-up quark polarization as to the gluon 
polarization. Including a degradation of the signal from finite (70%) beam polarization, a 
0.5% measurement of A~[ L would be required to determine Puix^) with an error of about 
5P u (x 2 ) ~ ±0.1. 

It is interesting to consider the complimentary spin asymmetry with transversely po- 
larized beams, Ajvtv- The formulas differ from that of the longitudinal asymmetry in two 
ways. First, the observed transverse asymmetry is also proportional to sin(</?), with ip the 
azimuthal angle of the production plane relative to the transverse polarization axis. Second, 
there is no contribution to the transverse asymmetry from processes with gluons, there- 
fore the leading source of transverse asymmetry should be the anti-up quark annihilation 
process, which will have a parton level analyzing power of —1. In the case of transverse 
polarization, none of the quark polarizations are known so detailed predictions are not very 
constrained, however, the sensitivity to polarization could be similar to what is described 
above for longitudinal u polarization in this region. 

VI. ACKNOWLEDGMENTS 

This work has been supported in part by grant RFFI # 95-02-05061 from the Russian 
Academy of Science and the National Science Science Foundation. 



11 



REFERENCES 



[1] HI Collab., S. Aid et al, Nucl. Phys. D470, 3 (1996). 

[2] ZEUS Collab., M. Derrick et al, Z. Phys. C69, 607 (1996); preprint DESY 96 - 076 

(1996), Z. Phys., to be published. 
[3] D. W. Duke and J. F. Owens, Phys. Rev. D30, 49 (1984). 
[4] H. L. Lai et al, Phys. Rev. D51, 4763 (1995). 
[5] J. Botts et al, Phys. Lett. B304, 159 (1993). 
[6] J. G. Morfin and W. K. Tung, Z. Phys. C52, 13 (1991). 

[7] A. D. Martin, R. G. Roberts and W. J. Stirling, Phys. Lett. B387, 419 (1996); ibid. 
B356, 89 (1995); ibid. B354, 155 (1995); Phys. Rev. D50, 6734 (1994). 

[8] F. Halzen and D. M. Scott, Phys. Rev. D18, 3378 (1978). 

[9] J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 
[10] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 
[11] CDF Collab., F. Abe et al, Phys. Rev. Lett. 71, 679 (1993). 

[12] CDF Collab., F. Abe et al, Phys. Rev. Lett. 73, 2662 (1996); ibid. 68, 2734 (1992). 
[13] CDF Collab., L. Nodulman, 28th International Conference on High Energy Physics 

(ICHEP '96), Warsaw, Poland, July 1996; FERMILAB - Conf - 96/337 - E (1996). 
[14] J. Huston et al, preprint MSU - HEP - 41027, CTEQ - 407 (1995). 
[15] C. Bourrely et al, Phys. Rep. 177, 319 (1989). 

[16] SLACE- 80, M. J. Alguard et al, Phys. Rev. Lett. 37, 1261 (1976); ibid. 41, 70 (1978). 

SLAC E - 130, G. Baum et al, Phys. Rev. Lett. 51, 1135 (1983). 
[17] EMC, J. Ashman et al, Phys. Lett. B206, 364 (1988); Nucl. Phys. B328, 1 (1989). 



12 



TABLES 



Ax 2 


.00- .05 


.05- .10 


.10- .15 


.15- .20 


.20 - .25 


.25- .30 


.30- .35 


.35 - .40 


.00 - .05 


426 


1,640 


1,460 


1,060 


660 


355 


173 


85 


.05 - .10 


1,260 


812 


264 


112 


51 


26 


13 


7 


.10 - .15 


365 


154 


45 


20.0 


9.0 


5.2 


2.9 


1.8 


.15 - .20 


70 


42 


13.0 


4.8 


2.5 


1.5 


0.64 


0.49 



TABLE I. Pythia cross sections in picobarns of the direct photon and jet production process. 
The statistical errors for the entries in these tables correspond to 35 events per picobarn. 
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44 


22 
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.05 - .10 


124 


116 


44 


20 


8.4 


4.7 


2.7 


1.2 


.10 - .15 


37.64 


23.07 


7.36 


3.26 


1.89 


1.25 


0.73 


0.20 


.15 - .20 


7.04 


6.46 


1.92 


0.73 


0.26 


0.23 


0.20 


0.15 



TABLE II. Cross sections in picobarns from the annihilation subprocess. 
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0.64 


0.35 



TABLE III. Cross sections in picobarns from the Compton scattering subprocess. 
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x\ and X2 


a+ 
(pb) 


a 

(pb) 


e 


C2 


Cl 


Co 


d\ 


d 


A(xi, x 2 ) 


6A/A (%) 
PYTHIA 


6A/A (%) 
STAR run 


0.175,0.075 


18.7 


13.0 


0.13 


0.68 


0.59 


0.27 


0.34 


0.15 


1.44 


3.8 


1.3 


0.225,0.075 


9.0 


5.6 


0.03 


0.75 


0.49 


0.29 


0.28 


0.13 


1.61 


5.5 


1.9 


0.175,0.025 


33.9 


19.0 


0.12 


0.70 


0.49 


0.31 


0.20 


0.13 


1.78 


2.9 


1.0 


0.225,0.025 


10.4 


5.5 


0.06 


0.84 


0.44 


0.33 


0.19 


0.13 


1.87 


5.4 


1.9 



TABLE IV. Cross sections and prediced asymmetry A(x±,X2) parameters as defined in Equa- 
tion The parameters are calculated with subprocess cross sections from Pythia and with d<i = 1 . 
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FIGURES 
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FIG. 1. The cos(i? 7 ) distributions at x\ = 0.175 and X2 = 0.075 from the reaction pp — > 7 + jet 
at center of mass energy 500 GeF and pt > 10 GeF in the acceptances of — 1 < t/ 7 < 2 and 
— 1 < 77jet < 2 for photon and jet respectively. Solid curve - Compton scattering subprocess, 
dashed curve - pair annihilation subprocess. The numbers in the vertical axis are cross sections in 
picobarns. 
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FIG. 2. Events distribution on a (r/ 7 x 77^) plot from the reaction pp — > 7 + jet at center of 
mass energy 500 GeV and ^ > 10 GeV in the acceptances of — 1 < t/ 7 < 2 and — 1 < 17^ < 2 for 
photon and jet respectively. The events are from x\ = 0.175 ± 0.025, %2 = 0.075 ± 0.025. In the 
top part of the figure events are from cos(# 7 ) < —0.7 and in the bottom part of the figure from 
cos(tf 7 ) > 0.7. 
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FIG. 3. Momentum fraction x\ (top) and X2 (bottom) distributions from the reaction 
pp —> 7 + jet at pp center of mass energy 500 GeV and pt > 10 GeV in the acceptances of 
— 1 < ?7 7 < 2 and — 1 < rjj e t < 2 for photon and jet respectively. The numbers in the vertical 
axes are cross sections in picobarms. Solid curves - total contribution from pair annihilation and 
Compton scattering subprocesses, dashed curves - contributions from the (qq) pair annihilation 
subprocess, and dotted curves - contributions from the Compton scattering subprocess. 
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FIG. 4. Two dimentional {x\ x X2) plots for events distribution from the reaction pp — ► 7 + jet 
at center of mass energy 500 GeV and pt > 10 GeV. The acceptances of the photon and jet are 
— 1 < r/y < 2 and — 1 < rjj e t < 2, respectively. The events angular distributions from the four 
dotted squares with the centers in the points (x 1 ,x 2 ) = (0.175,0.075), (0.225,0.075), (0.175,0.025) 
and (0.225, 0.025) are discussed in this paper. 
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FIG. 5. The r(x2) dependences of the forward - backward asymmetry A(x\, X2\ r(x2)) at 
x\ = 0.175 and X2 = 0.075 from the reaction pp — > j+jet at pp center of mass energy 500 GeV and 
Pt > 10 GeV in the acceptances of — 1 < ry 7 < 2 and — 1 < r/j e t < 2 for photon and jet respectively. 
The number under the dashed curve is the value of the forward-backward asymmetry at the limit 
r(x 2 ) > 1. 
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FIG. 6. The distributions of u, u quarks and G gluons in the proton multiplied by x (top) 
and functions r(x) and R(x) at Q 2 = 100 GeV 2 from CTEQ2L parametrization (bottom). In 
the top figure: solid curve - xu(x) quark distribution, dotted curve - xu(x) antiquark distribu- 
tion, dash-dotted curve - xu v (x) valence quark distribution and dashed curve - xG(x) gluons 
distribution; the the bottom figure: solid curve - r(x) and dashed curve - R(x). 
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FIG. 7. The r(x2) dependences of the relative error 8r(x2)/r(x2) in percentages at different 
points of x\ and xi'- top - x\ = 0.175, xi = 0.075; bottom - x± = 0.175, X2 = 0.025. Solid curves - 
relative error 5r{x2) /r{x2) related with uncertainty of R(x), if SR(x) / 'R(x) = 10%, dashed curves - 
uncertainty on r(x2) related with statistical error from asymmetry. The value of the relative error 
of forward - backward asymmetry 5A/A = 1%, which would correspond to a standard STAR run. 
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